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Abstract 

In unification models based on SU(15) or SU(16), baryon number is part of the gauge symme- 
try, broken spontaneously. In such models, we discuss various scenarios of important baryon 
number violating processes like proton decay and neutron-antineutron oscillation. Our analysis 
depends on the effective operator method, and covers many variations of symmetry breaking, 
including different intermediate groups and different Higgs boson content. We discuss pro- 
cesses mediated by gauge bosons and Higgs bosons parallely. We show how accidental global 
or discrete symmetries present in the full gauge invariant Lagrangian restrict baryon number 
violating processes in these models. In all cases, we find that baryon number violating inter- 
actions are sufficiently suppressed to allow grand unification at energies much lower than the 
usual 10 16 GeV. 



1 Introduction 



In all gauge theories reasonably verified by experiments, fermions transform as the fundamental represen- 
tations of the non-abelian gauge groups. Quarks transform as the fundamental representation of the color 
group SU(3), left-handed fermions are fundamental representations of the electroweak SU(2). Inspired by 
this, it is intriguing to consider the idea that all fermions transform like the fundamental representation of 
the grand unified gauge group. This leads to grand unified models based on the maximal symmetry group 
|Q] for each generation, SU(16), where the fermions all appear in the fundamental multiplet: 

^ L = (u r UbU y d r df,d y d r dbd y u r UbUy v e e~ e + u e ^j . (1-1) 

The indices r, b, y are three colors, and hats denote antiparticles for any fermion field ip: 

£=C 7 V, (1.2) 

where Cj^C^ 1 = — jT. Thus, for example, V e i is the antiparticle of the right handed neutrino v e R, 
assuming that it exists. The same pattern is repeated for other generations. Mirror fermions are needed 
to cancel the anomalies. One important feature of this model is that both baryon number (B) and lepton 
number (L), which are known symmetries of low energy physics, appear as gauge symmetries at high 
energy. In fact, this was one of the main motivations of Pati, Salam and Strathdee who first introduced 
such models Q. 

A new variant of these models has received some attention lately, where the gauge group is SU(15) |2|>[3|]. 
The difference with SU(16) is that the right handed neutrinos, which are not confirmed experimentally, are 
assumed not to exist, so that there are only 15 left-chiral fermionic fields per generation. Baryon number 
is still part of the gauge symmetry although lepton number is not. 

The interesting point about these models is that their characteristics are very different from the standard 
unification models based on the gauge groups SU(5), SO(10), Eq etc. For example, it has been shown that 
renormalization group analysis of certain symmetry breaking chains of these models yield low unification 
scales [|, |, |, g, as low as 10 8 GeV in some cases. All the known chains with such low unification scale have 
the property that they all break the unified group in such a way that at intermediate scales, quarks and 
leptons transform under separate subgroups of the gauge group. Because of low unification scales, these 
models do not suffer from the cosmological monopole problem ||, in sharp difference with SU(5) models. 
Important and interesting constraints on rare processes can be put in these models [f|, 0]. Although many 
of these points were first made with the SU(15) gauge group, it is now known that there are symmetry 
breaking chains of the original SU(16) gauge group as well which show these characteristics [||, p|, 10 [ . 



One crucial question arises now. How can a low unification scale be consistent with known bounds 
on proton lifetime? Of course, it is easy to see that gauge interactions do not violate B in the unbroken 



phase. This is another important difference with SU(5), SO(10) or E§ models. With a limited number of 
Higgs multiplets, it was argued that baryon number symmetry (B) is not violated || even after symmetry 
breaking has taken place. Subsequently, it was emphasized [||, |ll]] that since B is part of the gauge 
symmetry of these models, it must be broken spontaneously in order to avoid a massless gauge boson 
corresponding to an unbroken B symmetry, and therefore the Higgs sector must be expanded. 

Once this was pointed out, various scenarios of proton decay were considered by different authors 
[12, ||, ||]. Particularly powerful is the method of effective operators, which will be explained in Sec. [2|. 



A simple dimensional analysis performed with these effective operators [13, || shows that proton decay 
amplitude in these models are suppressed by as many as the fifth power of the grand unification mass, 
as opposed to the second power in the case of standard unification models. Because of this reason, low 
unification scale is consistent with the known bounds on the proton lifetime. 

However, this analysis was performed only with proton decay mediated by Higgs bosons, in a handful 
of scenarios. But there is another kind of contribution to the proton decay amplitude in these models. It 
is true that in the unbroken theory, each gauge boson carries a well-defined baryon number and therefore 
cannot mediate l?-violating processes |2|, [|. However, once the gauge group is spontaneously broken, 
gauge bosons with different baryon numbers can mix with one another and therefore the mass eigenstates 
of gauge bosons are not, in general, eigenstates of baryon number. They can therefore mediate baryon 
number violating processes. Although such contributions were discussed in some detail 0, m some 
early papers on the SU(16) model, the possibility of low unification scale was not realized at that time. 
In the context of low energy unification, it was discussed briefly only at the tree level in a very specific 
scenario ||. 

In the present paper, our focus is threefold. First, we discuss not only proton decay which is a \SB\ = 1 
process, but also neutron-antineutron oscillation which is a \SB\ = 2 process to see whether both are 
consistent with low energy unification. Of course, even when we will be explicitly talking about "proton 
decay", the comments can easily be translated to the baryon number violating decays of the neutron as 
well. Second, we include both gauge boson mediated and Higgs boson mediated processes in our analysis. 
Third, we point out accidental global symmetries in the full gauge invariant Lagrangians of the models 
which seriously constrain possible baryon number violating processes. For example, in one case we find 
that there are no operators involving four fermions which can give baryon number violation even after 
baryon number symmetry is spontaneously violated. In such cases, we extend the analysis of proton decay 
to operators involving six fermions, which has not been done before. 
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2 General considerations 



From the symmetries of the standard model alone, one can argue that the dimension of any proton decay 



operators must be six or higher 17], whereas for \AB\ = 2 processes it is at least nine [18]. Thus, we 
should be looking at non-renormalizable operators generated by the theory. In general, these operators can 
involve both ordinary and mirror fermions since the physical up quark, for example, can be a superposition 
of the ordinary and the mirror quark fields. However, for the sake of simplicity, we will assume that 
the mirrors do not mix with ordinary fermions. This can be attained naturally if we impose a discrete 
symmetry, 

xZ,(M) _> _\jf(M) wnere y(M) stanc i s for the mirror fields. This is the most popular discrete 
symmetry considered in the context of mirror fermions. In the presence of this symmetry the model ceases 
to remain vectorial and hence the fermion masses are protected. Hence the survival hypothesis is applicable 
to these theories. The mirrors can be heavier than the ordinary fermions if their Yukawa couplings are 
consistently larger. In that case, they will not figure in any of the low energy processes we will be discussing. 
We will also assume that there are no hitherto unknown bosons lighter than the nucleon mass. Low energy 
operators involving nucleons should then involve ordinary fermionic fields only. 



To analyze these operators, we adopt the procedure used in Refs. [^, 13 1, where one first constructs 
effective operators which are invariant under the gauge group. For baryon number violating processes at 
low energy, the full gauge invariant operators will in general also contain some scalar fields so that, when 
these scalars develop vacuum expectation values (VEVs), one obtains operators involving the fermionic 
fields only. If the VEVs are baryon number violating, the fermion field operator generated after putting the 
VEVs would violate baryon number. This is the main difference with operator analysis for proton decay 



performed in the context of SU(5) or SO(10) grand-unified models [16, 17], where the operators have to 
obey only the symmetries of the standard model since the unbroken grand unified model does not conserve 
baryon number. Here, the unbroken operators must obey the symmetries of the unified model, which of 
course is much larger than that of the standard model. This requirement severely restricts the type of 
baryon number violating operators that one can construct. 

The above discussion is applicable equally for any baryon number violating process induced by Higgs- 
boson and gauge-boson exchanges. We now consider proton decay in particular. Here, one needs an 



operator where the number of fermionic fields is at least four [16], and in most part we will discuss 



operators where the number is in fact four, except in Sec. 3.2 where we will find that such operators are 
forbidden for proton decay. 

For Higgs-boson mediated processes, the relevant bilinears are (^l) T C^l, where C is the conjuga- 
tion matrix for fermions. For gauge-boson mediated processes, the relevant bilinears are ^lIx^l- The 
important difference is that the first ^ appears with a complex conjugation here. If we put in the gauge 
indices, this will mean a lower index for the first and an upper index for the second one. Thus, the 
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gauge invariants constructed will have a quite different nature than the ones for the Higgs-boson mediated 
processes. 

The discussion so far implies that, for processes mediated by gauge bosons, the 4-fermionic operators 
must appear in the effective operator in the combination 



(2.1) 



where i,j,k,l denote gauge indices, and we have suppressed the generation indices. On the other hand, 
for Higgs boson mediated processes, the combinations should be 



(2.2) 



One can think of other operators like [($ i? ) i C($ fl ) i J R ) k C($ R )i\ or [(^ L yC(^ L ) j ] 
where is the multiplet which contains the antiparticles of the fields in ^l, but these are either just 
hermitian conjugates of the operators in Eqs. ( |2.1[ ) and (|2,2|), or can be Fierz transformed to them. So, we 
need not discuss them separately. 

Our goal is to find gauge invariant operators which can give rise to the four-fermion operators of Eqs. 
( |2.1| ) and ( |2.2| ) after symmetry breaking. This discussion involves the Higgs content of the model and the 
precise way in which baryon number is violated, and therefore has to be done separately for SU(15) and 
SU(16). This will be done in the ensuing sections. Also, as we said before, we will encounter specific 
models where 4-fermion operators are inconsistent with the symmetries of the model. For such cases, we 
describe here the general case where the fermionic part of the operator has 2n number of fields. Let us 
denote such an operator symbolically as K,f2 n \ij) 2n . The co-efficient JC(2 n ) has a mass dimension 4 — 3n. 
Thus, neglecting the masses of all decay products, a simple dimensional analysis will give 



mT' 7 }C 2 



(2n) 



Then, the known limits on the proton lifetime,^] 



r p > 3 x 10 32 yr 



implies 



IC {2n) <W~ 32 m*- 3n . 



(2.3) 



(2.4) 



(2.5) 



In the specific case where n = 2 (which is the most frequent case so that we will omit the subscript of JC 
in this case), we obtain 



/C<1(T 32 GeV~ 2 . 



(2.6) 



1 The precise bounds depend on the specific decay mode. For rough estimates, we use the same bound for all modes. 
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In standard unification models like SU(5) or SO(10), K, ~ g 2 Mg 2 , so that one needs Mc/g > 10 16 GeV in 
order for the models to be phenomenologically viable. In the models that we consider, we will see that /C 
is further suppressed by ratios of different mass scales, and that is why smaller unification scales will be 
consistent with phenomenology. 



3 Scenarios of baryon number violation in SU(15) models 



Various scenarios of baryon number violation has been discussed in the literature |l2j, |13|, in the context 
of the SU(15) gauge group. We present some such chains later in this section. In all these chains, for 
all symmery breakings above the weak scale, we use Higgs bosons either in completely antisymmetric 
representations, or in the adjoint representation or representations which can be obtained by taking tensor 
products of two adjoint representations. This is done for the sake of economy and definiteness. We also 
assume that, unless mentioned otherwise, the only Higgs boson multiplets present in the model are the 
ones which have VEVs. 

At the weak scale, however, we make an exception. Here, unless otherwise specified, we assume that the 
symmetry breaking is performed by a symmetric rank-2 multiplet S. This is motivated phenomenologically. 
If we use antisymmetric tensor to be the only field to couple to fermions, the fermion mass matrices would 
be antisymmetric. For three generations, this will imply that one mass eigenvalue is zero and the other 
two equal for particles of any given charge. This is very unrealistic, so we will not consider this possibility 
further. 



3.1 Baryon number violated by an antisymmetric rank-3 multiplet 
3.1.1 Symmetries of the model 



For the SU(15) model, Pal [13] introduced the most economic Higgs boson spectrum that leads to break- 
ing pattern with "un-unified" intermediate stages. The multiplets necessary for this purpose are the 

antisymmetric rank-3 multiplet ^ % i k \ the adjoint T l j, the symmetric rank-2 multiplet which gives 

Hi] 

fermion masses, and an additional one, H^i, which will be called the antisymmetric bi-adjoint since it 
appears in the tensor product of two adjoint representations and both the upper and the lower indices 
are antisymmetrized. Here and henceforth, the square and curly brackets denote antisymmetrization and 
symmetrization of indices. 

In Fig. [j], we show the complete chain of symmetry breaking, where the numbers n denote a factor 
SU(n) in the gauge group if n > 1, and a U(l) factor if n = 1. Thus, at the highest stage, the multiplet 
§[ijk] develops a VEV <V ee ~ e+ J>, which breaks the unification group SU(15) down to SU(12) 9 x SU(3) € , 
where the subscripts q or i denote that only quarks/antiquarks or leptons/antileptons are non-singlets 
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under the subgroup. This VEV also breaks lepton number by 1 unit. At the next stage, SU(12) breaks 
to SU(6) 9L x SU(6) gJ j x U(1) B . This can be performed by a VEV in the adjoint representation, as shown 
in the figure and explained in the figure caption. At the scale M^l, the SU(6) g L breaks to its maximal 
subgroup SU(3) ® SU(2), under which the fundamental of SU(6) transforms like (3,2). In SU(6), the lowest 
dimensional multiplet which has the component whose VEV can induce this breaking is the 189-dimensional 
antisymmetric bi-adjoint. Naturally, it is contained in the antisymmetric bi-adjoint representation of the 
SU(15). At the next stage, SU(6) gJ? breaks to SU(3) uJ? x SU(3) di? x V(l) qRA . The SU(3) factors here 
operate non-trivially only on the ul and di components respectively, whereas the U(l) quantum numbers 
are defined to be +1 for the up- type quarks and —1 for the down- type ones. Baryon number is broken at 
the next stage, where SU(3) ufi x SU(3) di? also is broken to the diagonal SU(3) subgroup which is called 
SU(3) 9 fl. At this stage, the gauge group that appears is the square of the standard model gauge group, 
where the quarks and leptons transform under different SU(2) and different U(l) factors. This has been 
discussed under the name "un-unified" model by some authors p0[| . On the other hand, the right and left 
chirality of quarks transform under different color groups, which has been discussed in the literature under 



the name "chiral color" [15, 21]. At the next stage, the standard model gauge group appears, which is why 
this scale is called Ms- 

It should be understood that some variations of this chain are obviously possible. For example, the 
scale M^ q L can be lower than M^r or even M^. On the other hand, some scales can merge, so that the 
standard model is reached in less number of steps. These will not essentially change the conclusions of 
the subsequent discussions and hence will not be discussed separately. Similar comments apply for other 
chains which we will discuss later in the paper. 

For this and various other scenarios that we are going to discuss, we find that the full gauge invariant 
Lagrangian involving the specified fields often contains some accidental global or discrete symmetries 
which commute with the gauge symmetry. These restrict the type of potentially baryon number violating 
operators. To see such symmetries in the present case, let us write the full Lagrangian in the following 
suggestive manner: 

C = C Q + C. (3.1) 

Here, Co is the part which is invariant under independent phase rotations of all complex multiplets present 
in the model. These would include, e.g., all gauge interactions, scalar interactions involving only the 
adjoint Higgs multiplet and the antisymmetric bi-adjoint. There will also be some terms involving other 
multiplets, e.g., terms like ^ k ^ij^ or S lJ SjkS kl Su. Obviously, symmetry of Cq is much larger than the 
gauge symmetry. However, £' contains other terms which are allowed by the gauge symmetry. In the 
present case, the Yukawa couplings are the only terms which fall in this class. Thus, 

C! = y '(* L ) k (*L) 1 Ski + h.c, (3.2) 
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where the generational indices on $i and y have been omitted. However, it is easy to see that, even 
with this term, the full gauge invariant Lagrangian has the following accidental global charges which are 
conserved: 



Multiplet 




s kl 




Qi 


1 


2 





Q2 








1 



(3.3) 



Notice that the global phase of the multiplet $ is indeed a global symmetry of the Lagrangian [19]. In 
addition, there is another one, which has been labeled as Q\. 
Consider now a generic effective operator of the form 

which stands for 2f number of fermionic fields with upper indices, ns number of the multiplet S with 
upper indices etc. Each of these numbers can be positive, negative (if the relevant multiplet contributes a 
net number of lower indices) or zero. Thus, for example, the operator $* J '*<i>yk will have = since 
contributes n$ = 1 but contributes n$ = —1. On the other hand, $^ fc $ Zmn SuSj m Sk n has n$ = 2, 
ns = -3. 

Conservation of the charges Q\ and Q2 tells us that, in Eq. (|3T 



f + n s = 0, (3.5) 
n$ = 0. (3.6) 

Any effective operator generated by the theory must then obey these two conditions, and we will discuss 
some such operators below. Note that both these conditions remain unaffected if the operator in Eq. (^ 



It il 

contains the adjoint or the multiplet , since they contribute equally to the number of upper and lower 



indices, and since they are neutral under the global symmetries of Eq. ( |3.3| ). 

One general characteristic of baryon number violation in this model can be immediately noted. As we 
said earlier, in order to obtain purely fermionic operators, we need to replace the scalar fields in Eq. ( |3.4|) 
by their VEVs. From Fig. |l], we note that $ has three types of VEVs. One of these gives 5B = —1, and 
each of the other two give 5L = 1. Let us say that in the operator with fermionic fields only, there are 
n$ B VEVs of the first type, and n$ L VEVs of the second kind. Since the purely fermionic operator comes 
from the gauge invariant operator of Eq. (3.4), they must obey Eq. (3.6), which implies n<j> s + n$ L = 0. 



The total violation of B and L in the purely fermionic operator can now be written as 

5B = — n$ s , 5L = n$ L = — n<j> s , (3-7) 

so that 

6(B-L) = 0. (3.8) 



This immediately tells us that there is no n-n oscillations in the model. 

One comment needs to be made here. Lepton number is not part of the gauge symmetry in SU(15). 
However, it is well-defined for all components of ^f. This can be used to assign lepton number to the gauge 
bosons and Higgs bosons. In other words, for any multiplet we can count a lepton number +1 for 

each occurence of the indices 13 or 14, and —1 for each occurence for the index 15. Since lower indices 
are complex conjugates, they will have just the opposite assignments. It is then easy to see that lepton 
number conservation is assured in the Lagrangian by gauge invariance, and therefore must be violated 
spontaneously. This is one characteristics of this particular version of SU(15) models which is not shared 
if we break baryon number by higher rank multiplets, as we will see later. 

3.1.2 Proton decay operators 

Since proton decay requires both baryon number and lepton number violation, and since both these vio- 
lations come through VEVs of different components of we must need at least two factors of <5 in the 
gauge invariant operator. Of course, one of them must come with upper indices and the other with lower 
indices in order that Eq. ( |3.6D is satisfied. 

Gauge boson mediated : For gauge boson mediated proton decay for which / = in the generic 
operator of Eq. (|3.4| ), we obtain n$ = from Eq. Indeed, there is one operator which is consistent 

with all these numbers: 



(*L)a A (* L y (*L)k7x(*L) 1 ^^jlr, (3.9) 



In Fig. ^, we have shown a tree-level diagram which gives rise to a particular component of this operator. 
This is the component with {ikr} = {udd}. Since the other VEV must have one index contracted with 
this one, and it has to be lepton number violating, we have to use either {jlr} = {due} or {ddv e }. These 
two possibilities are shown in Fig. |2|a and respectively. Once the VEVs are put in, 4-fermion operators 
result, whose co-efficient can be easily computed by looking at these diagrams: 

g 2 M B M s f . 

/Cl " mP^ • (3 - 10) 

Here, the factor g 2 is just the gauge coupling constant coming from two vertices with fermions. Mgjg 
and Ms/g give the VEVs, and the four-boson vertex gives a factor g 2 . The denomintor comes from the 
propagators of the gauge bosons. In both diagrams, the gauge boson coming out of the left vertex has the 
quantum numbers of a diquark. Such gauge bosons belong to the coset space SU(12) (J /[SU(6) (?i x SU(6) ?i? ], 
and acquire masses of order of SU(12) 9 breaking scale, Myi- The other gauge boson which couples a quark 
to a lepton belongs to the coset space SU(15)/SU(12) (J x SU(3)£, and therefore has mass at the unification 
scale Mq. 
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Proton lifetime bounds now imply, from Eq. (|2.6D, the constraints 



Ml ■ M * 2 > 10 32 GeV 2 . 



(3.11) 



g 2 M B M s 

Notice that Mb and M$ are by definition smaller than M12, and g < 1. Thus this condition can be satisfied 
with a low grand unification scale. 

We now discuss the proton decay modes obtained from the diagrams of Fig. ||. In the figure, we have 
suppressed all generation indices of the fermions. We now notice that there is a property of the operator 
of Eq. Q3.9|) which forbids all fermions to be of the same generation. This is because 



(*Z)i(*fl)fc] [i^) l (^L) j 

2 ($ R ) i c($ fl )J U* L ) l c(* L 



(3.12) 



Here, the first step is obtained by the definition of and the next one is obtained by Fierz transforma- 
tion. Since the matrix C is antisymmetric, the last form shows that the spinor indices of W and ^ k are 
antisymmetric in Eq. (3.9). Because the gauge indices are contracted with 3> lfcr which is antisymmetric, 
the gauge indices i and k are also antisymmetric. Therefore, in order to satisfy the Fermi principle, they 
must be antisymmetric in their (unshown) generation indices. The same comment can be made about ^ 
and Vl/'. Thus, disregarding charm quark fields which will be kinematically forbidden in proton decay, we 
obtain that the 4-fermion operator generated by Fig. |2|a is 



(3.13) 



Similarly, Fig. ^b generates 



(3.14) 



The first one predicts the decay mode 



(3.15) 



whereas the second one gives 



P 



(3.16) 



Notice that these are unusual decay modes which are suppressed in unification models like SU(5) or SO(10), 
although they occur in their supersymmetric versions. In the present case, it is predicted in absence of 
supersymmetry because of the Fermi symmetry between all particles in a generation. 
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Higgs boson mediated : In this case, we should put / = 2 in Eq. ( |3 . 4| ) , as discussed in Sec. [2|. The 
constraint of Eq. (|3.5| ) now implies n$ = —2. Operators of this type were discussed earlier by one of us 
[13]. Here is one example: 



2 = [^LyC(^ L y] [^ L ) k C^ L ) l ]<S> ikr ^ r S lp S jq . (3.17) 

Since the effective operator now involves two occurrences of the field S whose VEVs are of order Myy, any 
contribution coming from this operator must have a suppression factor of (Mw/Mq) 2 . For example, the 
tree diagram of Fig. || gives a contribution 

K2 ^4«, ( , 18) 

Here, the quantity is the mass of a typical fermion, and comes from the Yukawa couplings. The quartic 
scalar couplings are denoted by A55 and As$ in an obvious notation, and we have assumed that all the 
virtual colored scalars in this diagram have masses of order Mq, the largest scale in the model. Then 
notice that 

Even if Ass,As$ ~ 1 which is the limit allowed by perturbative procedure, we obtain JC2 <C K-i since 
M\2 < Mq by definition and raj Mq. 

3.2 Baryon number violated by an antisymmetric rank-4 multiplet 
3.2.1 Symmetries of the model 

The rank-4 antisymmetric representation A was introduced by Brahmachari, Sarkar, Mann and Steele 
(BSMS) ||. A possible chain of symmetry breaking has been shown in Fig. f|. Notice that the symme- 
try breaking above the scale Mb is performed by the same VEVs as in Fig. Q. At the scale Mb, the 
still unbroken gauge group ^qL^qL^uR^dR^-B^qRh^iL^-tY is broken by the VEV of the rank-4 multiplet, 

A uude+ _ The and the 3dR gub 

groups combine to give the diagonal subgroup which we call 3 q R. Also, 
out of the three U(l) factors, one combination breaks, leaving two unbroken ones, one of which is the 
hypercharge of the standard model, and the other is called XJ(1)f, whose generator is proportional to 
diag (9 {6) , -22(3), 4 (3) , -7 (2) , 14). 

Notice that the VEV A. uude breaks both baryon number and lepton number by —1. Since this is the 
only baryon number violating VEV in this scheme, it must appear in the gauge invariant operator giving 
rise to proton decay. Moreover, it must appear an odd number of times. Thus, it is obvious that if the 
model contains a discrete symmetry 

A -> -A (3.20) 
11 



with all other fields invariant, one cannot generate any term that violates baryon number by an odd 
integer. Proton will then be absolutely stable. This comment applies irrespective of whether proton decay 
is mediated by gauge boson or Higgs boson exchange. 

Even if such a symmetry is not imposed on the Lagrangian, analysis of the full gauge invariant La- 
grangian reveals accidental global symmetries as discussed in Sec. 3.1. H restricting the type of potentially 



baryon number violating operators. To see this, we use the notation of Eq. ( |3.1| ) and note that here 

£' = y(y L ) k (V L ) l S kl + XA klmn S pr <^ M ^ mnr + A'[AAA$] e + h.c. , (3.21) 

where in the last term, the indices (not shown) are all upper indices which are contracted by an antisym- 
metric e-tensor having 15 indices, which is indicated by the square bracket with a subscript e. The first 
term is the Yukawa coupling term. 

It is easy to see that, even with the presence of the above terms, there is a global U(l) symmetry of 
the Lagrangian under which the quantum numbers of various multiplets are as follows: 



Multiplet 



Charge 



(3.22) 



1 9 — — — 

Consider now a generic effective operator of the form 

(f) 2/ ,S ns $ n *A nA , (3.23) 

in the notation used in Eq. ( |3.4[ ). The global symmetry of Eq. ( |3.22| ) implies 

2f + 2n s + -n$ - -n A = . (3.24) 

On the other hand, all the indices should be contracted, which means that the total number of upper 
indices should either be zero or be divisible by 15 (so that they can be contracted by e-symbols). For 
the model of Sec. [O], this condition is already contained in Eqs. Q3.5| ) and Q3.6| ). Here, it produces an 
independent condition^ 

2/ + 2n s + 3ncj, + 4n A = 15iV , (3.25) 

where iV is an integer, denoting the number of times a vertex involving the e-symbol appears in the diagram 
giving rise to the operator of Eq. ( |3.23| ). As noted in Sec. |3,l.l| , both conditions remain unaffected if the 
operator in Eq. ( |3.23 ) contains the adjoint or the antisymmetric bi-adjoint. 



The solution of Eqs. (|3.24|) and (|3~25| ) can be written as 



/ + n s = n A - 3N, n$ = 7N - 2n A . (3.26) 



2 This condition is necessary, but not sufficient, since it does not take into account the fact that the indices to be contracted 
by the e-symbols have to be antisymmetric. 
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Let us now check what the above solution means for the violation of baryon and lepton numbers. Baryon 
number, as noted before, is part of the gauge symmetry and is broken only spontaneously through the 
VEV of A. Thus, clearly, 

SB = -n A . (3.27) 
On the other hand, lepton number violation comes from three different sources: 

• each VEV of A (with upper indices) gives SL = — 1; 

• each VEV of (with upper indices) induces 5L = 1; 

• each occurence of a term with an e-symbol will have 15 upper indices which are all different, con- 
tributing to an explicit violation 5L = 1 in the unbroken Lagrangian. 

Taking all these contributions, we can write 

8L = -n A + n$ + iV . (3.28) 

Using Eqs. ( |3 . 2 C - 3728| ) , we therefore finally obtain 

5{W - L) = -8N , (3.29) 

which is the selection rule for this model. Immediately, it tells us that in this model, there cannot be any 
neutron- antineutron oscillations. 

3.2.2 Proton decay operators 

Specializing to the simplest case when N = 0, Eq. ( |3.29| ) tells us that 3-B — L is conserved, which means 
that there will be three leptons in the final state for proton decay. This cannot occur with four fermionic 
fields only, since three of these fields must be quark/antiquark fields in order to obtain a 5B = — 1 operator. 
For other values of N, one needs even higher number of leptons/antileptons in the final state, which cannot 
be accommodated in a 4-fermion operator for the same reason. Thus, we conclude that in this model, there 
is no proton decay operator with four fermionic field operators. The result is true for operators mediated 
by gauge or Higgs bosons. 

The lowest dimensional operators will thus have six fermionic fields. They can have / = 3 where all 
the indices are upper. Alternatively, they may have / = 1 where two of the fields have lower indices, but 
the other four have upper ones. Of course, one can similarly have / = —3 and / = —1. 



Among these possibilities, / = 1 can yield a solution to Eq. ( 3.26| ) with smallest number of scalar fields, 
given by 

/ = l, n s = 0, ra$ = -2, n A = 1. (3.30) 
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An operator of this type is: 



O 



(*i)i7A(*i) J (*£)*tW [(^LfC^LY) A ikab <S> ajr <S> plb . (3.31) 



We show in Fig. [| how this operator can arise at the tree level. The amplitude of the purely fermionic 
operators can be easily determined. Assuming the scalar interaction couplings to be of order unity, we 
obtain 

gM G M B M s ( . 

~ m^mJmI — mi ' (3 - 32) 

12 G due-e+ eu 

where the last two factors in the denominator represent the masses of the internal Higgs boson lines. Of 
these, the former one is a colored boson, whose mass is expected to be of order Mq. But the latter one 
is uncolored, whose mass we keep as an unknown. Experimental bounds, however, tell us that, being a 
charged scalar, its mass cannot be much less than 100 GeV. Thus, if this operator contributes to proton 
decay, using Eq. (2.5), we can rewrite Eq. (|3.32| ) as 



Since Myi > Mb, Ms by definition and g < 1, this bound can be satisfied for any unification scale larger 
than about 10 9 GeV. 

However, there is a subtle reason why this operator cannot contribute to proton decay. In order to 
accommodate baryon number violation, the indices on the fields A must be uude + in any permutation. 
Now, these indices contract either with the indices of <£, or those of ^f. But $ does not have any VEV 
which contains the index u. Thus, both the indices on Vl/ have to be u indices. However, as argued in 
connection with Eq. ( |3.12 ), the fields (&L,)i and (^L,)k must come from different generations. Therefore, 



one of them must be the charm quark and therefore proton decay is kinematically forbidden from this 
operator. 

To get out of this impasse, one can use a slightly modified operator: 



O' 



(*z)m(*x) 5 l(y L yC(* L y] T m A mkab <S>ajr$ P i b . (3.34) 



This is still an operator of the type of Eq. ( |3.30| ), but now the gauge indices i and k are not antisymmetric, 
and therefore and (^i)k can refer to fields from the same generation. A diagram for this operator 

can be obtained from Fig. [| by attaching an adjoint Higgs boson to any line which carries at least one 
SU(12)g index. This will provide further suppression to the 4-fermion operators since the extra propagator 
is expected to have a mass Mq, but the largest VEV available for the adjoint multiplet is at the scale M\%- 
The quark level transition induced by this operator is ue~e~v — > uu, which implies a decay mode 

p — > TT~e + e + u e . (3.35) 
14 



3.3 Baryon number violated by an antisymmetric rank-5 multiplet 



In an early paper, Frampton and Kephart 12] discussed baryon number violation by the VEV of an 
antisymmetric rank-5 multiplet jl l i klm \. Although less economical than the ones discussed above, we 
include this possibility for the sake of completeness. Fig. ^ gives a chain involving this rank-5 multiplet. 
The VEV that breaks U(1)b has the gauge transformation properties of dddv e e~ , i.e, it has B = —1, 
L = 2. Notice also that since this VEV does not involve both u and d type indices, it cannot break 

SuR&dR part of the symmetry to 3 q ji, as is done in the models described earlier. Therefore, this symmetry 

fill 

breaking is performed by a VEV in the antisymmetric bi-adjoint Hl^n- This multiplet certainly has a 
component which is the antisymmetric bi-adjoint of the subgroup Q q R and singlet under the rest. This part 
is a 189-dimensional representation of SU(6) which has a component that transforms like (8, 8, 0) under its 
subs roup ^uR^dR^qRK- A VEV here would perform the desired symmetry breaking. On the other hand, 
one now does not need the adjoint to break the 6 q R subgroup, since the baryon number violating VEV 
itself performs the job. In fact, the VEV I^J dddu ^ e ) also breaks the leptonic subgroup ?>t to 2?l, and the 
leptonic and quark hypercharges combine to the total hypercharge of the standard model. 
In this case, using the notation introduced earlier, we obtain 

£ = y^> L ) k ^> L ) l S k i + fi[JJJ] t + h.c. . (3.36) 

Obviously, the entire Lagrangian respects the discrete symmetry 

J e 2ni/3 J (3.37) 

with all other fields neutral. This is a Z3 symmetry. The number of J fields in any effective operator 
arising in this model must then be a multiple of 3. Since baryon number violation comes from the VEV of 
J only, in purely fermionic operators we will have 

\SB\ = 3N (3.38) 

for some integer N. Therefore, neither proton decay nor neutron- antineutron oscillation is possible in this 
modelJ^ Notice that this conclusion is reached only from the accidental symmetries present in the full 
gauge invariant Lagrangian. 

3.4 Introduction of antisymmetric Yukawa couplings 

So far, we have assumed that the only Higgs bosons which can couple to fermions belong to the symmetric 

rank-2 multiplet . The situation changes if, in addition there is also the multiplet which couples 
3 Indeed, the baryon number violating diagram given by Frampton and Kephart jli) for this model has SB = 3, as was 
noted by one of us earlier p^ ]. 
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antisymmetrically. In this case, some of the symmetries described in the above sections may be broken 
explicitly and hence more baryon number violating processes may be allowed. 

For the model of Sec. 3.1, such is not the case. We still have the condition in Eq. (3.6), which leads 



to B — L conservation. However, in the model of Sec. [3.2| , there is an important change. This is because, 
with the introduction of the multiplet A, there are the following new terms which are allowed in 

C' A = y A (*L) k (*L) l A kl + vA i3 A kl A ijkl + h.c. . (3.39) 



There is now no way that one can assign a quantum number of A which keeps the symmetry of Eq. ( 3.22j ). 
Thus, one can have the following operator: 

o = [(^ L yc^ L y] i(* L ) k c^ L ) l ]A lJkl . (3.40) 



In Fig. [5], we show how this can be generated through the interactions appearing in Eq. ( 3.3S| ). In the 



figure, we suppressed the generation indices. Turning to Eq. ( 3.40Q , we see that since the gauge group 
indices i and j appear in antisymmetric combination in A l ^ kl , and since the matrix C is antisymmetric, 
the generation indices for the two fermionic fields in the first bilinear must be antisymmetric in order to 
maintain Fermi symmetry. The same can be said about the fermionic fields in the other bilinear. Thus, 
the quark level operator coming from Fig. [7| is [2C/i + ] [uCs], which gives rise to a proton decay mode 

p -> fi + K° . (3.41) 

The amplitude for the 4-fermion operator is given by 

assuming, once again, that the colored Higgs bosons have masses of order Mq. The quantity y A in 
this formula stands symbolically for two factors of the Yukawa coupling with the multiplet A. Since the 
antisymmetric Yukawa couplings, if any, are expected to be smaller than the symmetric ones, and since 
Mb < Mq by definition, this again shows the suppression of proton decay rate. 

For the model of Sec. [T^, the changes are more dramatic. Here, the extra terms can appear in C due 
to the introduction of A are given by 

C! A = yA^L) k {^L) l A kl + \jj* A [JJ®A] t + fi'Jv klm <S> klm A l3 + h.c. . (3.43) 

There still is a Z3 symmetry in the full Lagrangian, defined as follows: 

1 

(3.44) 



Multiplet 


^jjfc gkl j^kl ^klm jijklm 


Z3 charge 


1 2 2 2 1 



Consider now a generic effective operator of the form 



(i$ffS n sA nA <f> niS 'J nj (3.45) 
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in the notation used before. Using the Z3 symmetry and the requirement that all indices must be contracted, 
we obtain the following conditions: 



2/ + 2{n s + n A ) + 2n$ + nj = 3N , 
2/ + 2{n s + n A ) + 3n$ + 5nj = 15JV' , 

where N and N' are both integers. Thus, 

n-j, = 15N' -3N - 4nj , 2/ + 2(n s + n A ) = 9N - 30N' + 7nj . 



Following arguments similar to those in Sec. 3.2, we now obtain 



so that, using Eq. ( p. 48 ), we obtain 



SB = -nj, 

5L = 2nj + U& + N' , 

5(2B — L) = 3N — 16N' , 



(3.46) 
(3.47) 

(3.48) 

(3.49) 
(3.50) 

(3.51) 



which is the selection rule in this case. 

For proton decay which requires nj = 1, notice that the integer iV must be odd because of its definition 
in Eq. ( |3.46 ). The solution of Eq. fl3,48|) involving minimum number of scalar fields is now given by 
N = -1, N' = 0, i.e., 



-1, / + n s + n A = -1 



(3.52) 



Using Eqs. ( 0.49[) and (^.50 ), it is now easy to see that in this case, proton decay operators will satisfy the 
selection rule 5(B + L) = 0. 

Example of a gauge boson mediated diagram of proton decay is provided in Fig. |8|, which has / = 
and n A = — 1. The operator here has the form 



(3.53) 



and Fig. || shows how it can be generated at the tree level. Because of the Fierz transformation property 
shown in Eq. ( |3.12j ), the generation indices of and \& fc must be different here. But the same cannot be 
said about and since their gauge indices are not antisymmetric. Thus, the quark-level transition 
obtained from Fig. S is de + — > sd. This implies a proton decay mode 



p — > ir + K + e 



(3.54) 



which conserves B + L, as argued before on general grounds. The coefficient of the 4-fermion operator is 
given by 

[i'M B M s M w 



£1 ~ 
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(3.55) 



assuming that the colored scalar internal line has a mass of order Mq. Once again, since Mb, Ms and 
Mw are each smaller than either M\2 or Mq by definition, a low unification scale is allowed. 

A Higgs boson mediated diagram, with / = — 2 and ha = 1, was given in Ref. [13]. The operator 
responsible for this is: 

2 =[{*L) i C(* L )>] [(* L ) fc C(* L ) z ] J ijkqr $™ r A lp . (3.56) 

It apparently looks like it has the same VEVs as the operator in Eq. ( |3.53| ). But this need not be the case, 
as seen from Fig. ||. Here, one can use the VEV of <3? which occurs at the scale Mq. Thus, we obtain for 
the strength of the 4-fermion operator 

\M W ) Ml V ' 

Depending on the magnitude of the scales Ms and M12, this may or may not dominate over the gauge 
boson mediated decay. The decay mode is the same as that given in Eq. ( |3.54 ) since Vl/ 1 and & have to 



belong to different generations in the operator of Eq. ( |3.56| ). 

It must also be noticed that, unlike the previous models, neutron-antineutron oscillations are not ruled 
out in this model. However, it is very suppressed. This can be seen from Eq. ( |3.51 ), where we can put 
5B = 2 and 5L = as is necessary for neutron-antineutron oscillations. The simplest solution for this 
situation is obtained when N = —4, N' = —1, which means 

nj = — 2, n<j> = 5, f + n-s + tia = — 10. (3.58) 

Obviously, it is a very high dimensional operator, so we will ignore it. 



4 Scenarios of baryon number violation in SU(16) models 

It might seem that SU(16) scenarios of baryon number violation should look similar to the SU(15) ones, 

since the groups are not all that different. There are, however, some important differences, which should 

carefully be taken into account. The first is that baryon number violation occurs spontaneously and 

therefore is sensitive to the choice of the Higgs sector, as amply demonstrated in Sec. [3|. Being a larger 

group, SU(16) in general requires more VEVs to break it down to the standard model gauge group, which 

affect the operator analysis. Secondly, baryon number processes like proton decay involves lepton number 

violation as well, and the latter is very different in the groups SU(15) and SU(16). The reason is that in 

SU(16) lepton number is part of the gauge symmetry and can be violated only spontaneously. This is a 

difference from the SU(15) models where the Lagrangian can violate lepton number. Thirdly, the symmetric 

rank-2 multiplet of SU(16), unlike its SU(15) version, can have a lepton number violating VEV S vu , which 

does not violate the symmetries of the standard model. This VEV can give neutrinos a Majorana mass at 

the tree level. To keep our discussion simple, we will neglect this VEV. 
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4.1 Baryon number violated by an antisymmetric rank-4 multiplet 
4.1.1 Symmetries of the model 

Breaking SU(15) down to SU(12) 9 xSU(3)^ requires the VEV of an antisymmetric rank-3 multiplet. Sim- 
ilarly, breaking SU(16) down to SU(12) 9 xSU(4)£ requires the VEV of an antisymmetric rank-4 multiplet 

A. Therefore, it is reasonable to try to see if there are suitable VEVs in the multiplet A, the adjoint T l j 

Hi] 

and the bi-adjoint H^J which can break the grand unification symmetry down to the symmetry of the 



standard model. In Fig. 10, we show how it can be done. At the weak scale, the symmetry is broken by 
the rank-2 symmetric multiplet , as before. 
For this model, we find 

£ = y^LTi^LySij + A'[AAAA] e + h.c. . (4.1) 

It is easy to see that it has an accidental global U(l)xi?4 symmetry, under which the charges of various 
multiplets are as follows: 



Multiplet 




s** 




U(l) charge 


1 


2 





Z4 charge 








1 



(4.2) 



Considering now a generic effective operator of the form 

($)2/5«S A nA i (43) 

The requirements of the U(l)xZ4 symmetry and of the contraction of all indices give the following condi- 
tions: 

2/ + 2n 5 = 0, (4.4) 

n A = 4N, (4.5) 

for some integer N. Notice that both baryon number and lepton number violation come from only one 
VEV, viz., (A"^). This VEV gives SB = -1, 5L = -1. Thus, in this model, 

S(B-L) = 0. (4.6) 

Once again, neutron- antineutron oscillation is not possible in this model. The possibilities of proton decay 
are discussed below. 

4.1.2 Proton decay operators 

Obviously, the simplest solution to Eqs. and flj.5|) are given by 

f = n s = n A = 0. (4.7) 
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An operator of this type is: 



(4.8) 



This can give gauge-boson mediated proton decay, as shown in Fig. [0]. Notice that this diagram is very 



similar to Fig. y. The analysis is also the same, leading to the constraint in Eq. ( 3.11 ). In fact, one can 
also show that Higgs boson mediated diagrams, having n$ = 2, will be suppressed in this model, as shown 
in Eq. 

4.2 Baryon number violated by an antisymmetric rank-3 multiplet 

Deshpande, Keith and Pal S advocated a model where baryon number symmetry is violated by a VEV 



of an antisymmetric rank-3 multiplet as in Sec. 3.1. In Fig. 12, we show this chain with some slight 
modifications which helps eliminate a fundamental Higgs multiplet which was used by them. 
With the introduction of the multiplet there is one more term in CJ : 



a = y^Ly^LYSij + AA Wmn s pr * Wi '* mnr + a'[aaaa] £ + h. c . . 

However, we now have an accidental U(l)xZs symmetry, with the following charge assignments: 



(4.9) 



Multiplet 




S# 






U(l) charge 


1 


2 


1 





Zg charge 








1 


2 



(4.10) 



So now, the generic effective operator of the form 

(*) 2/ S ns $ n *A" A 

is subject to the following constraints: 

2/ + 2n s + n$ = 0, 
n$ + 2n A = 8iV , 



(4.11) 



(4.12) 
(4.13) 



for some integer N. The first of these equations now implies that n$ must be even. The simplest so- 
lution to these conditions is given by all n's being zero, which gives the operator of Eq. (|3.9|), and the 



phenomenological conclusions are the same as in Sec. 3.1 



5 Conclusions 

We have analyzed a variety of symmetry breaking chains within the gauge groups SU(15) and SU(16) in 

which the grand unified gauge group breaks to "un-unified" subgroups under which quarks and leptons 
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have separate symmetries. As mentioned in the Introduction, such chains are interesting because some of 
them are known to predict low unification scales, sometimes as low as 10 8 GeV. Our analysis shows that 
low unification scale is not phenomenologically ruled out in these models because proton decay operators 
are very suppressed. The amount of suppression, of course, depends on the Higgs boson sector of the 
models and therefore varies from one model to another. We have also shown that in these models, since 
operator analysis can be performed on the full gauge invariant operators, and since such operators have 
a large Fermi symmetry, the proton cannot decay into non-strange hadrons. Such modes are preferred in 
supersymmetric SU(5) or SO(10) models, but here we obtain this conclusion without any supersymmetry in 
our models. In fact, inclusion of supersymmetry in SU(15) or SU(16) models typically make the unification 



scales large [22]. In that case, with all the suppression mentioned in this paper, proton decay should be 
unobservably slow. 

One remarkable result that comes out from our analysis is that, most of these models contain accidental 
global or discrete symmetries. Of course, these symmetries depend on the Higgs boson contents of the 
model, much like the B — L conservation in the simplest SU(5) unification model. Such symmetries 
provide selection rules to baryon number violating processes. For example, neutron-antineutron oscillation 
is strictly forbidden in most of the models, as we pointed out. Of course, one can always further complicate 
the models, using more Higgs boson multiplets than are necessary for breaking the symmetries. Presence 
of these multiplets will explicitly break some or all of the accidental symmetries that we discovered, and 
therefore will allow more baryon number violating processes. We provided examples of this by introducing, 



in Sec. 3.4, the antisymmetric rank-2 tensor which couples to fermions. Another example could be a 
multiplet X^ k ^ mn ^ , which exists in the symmetric part of the tensor product of two antisymmtric rank- 
3 multiplets. Once this multiplet is introduced, one can show that neutron-antineutron oscillations become 
allowed in most cases through the operator 

[(y L yc(v L y] [(* L ) k c(* L ) 1 ] l^L) m c(* L T}x {m[lmn]} . (5.1) 

Our analysis, however, deals mostly with "minimal" models in the sense that we do not introduce any 
Higgs boson multiplets which are not necessary for symmetry breaking. 

Our results can also be used to look for other baryon and lepton number violating processes in these 
models. For example, using Eqs. (p.8|), (|3.29|) and (|4.6|), we can conclude that neutrinos cannot have any 



Majorana mass in the SU(15) models of Sees. 3.1 and 3.2, as well as in the SU(16) model of Sec. 4.1. For 
the first of these models, this result was derived by earlier authors [19], but for the other models, the result 
is new. This and other new results can be readily derived from the accidental symmetries that we have 
discovered in this article for various models of interest. 
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Figure 1: SU(15) symmetry breaking where baryon number is broken by the VEV of an antisymmetric 
rank-3 multiplet. If one considers the adjoint Higgs multiplet T as a traceless matrix, its VEVs are diagonal 
and the notation 1( 6 ), e.g., stands for six consecutive entries of unity. In the multiplet <£, the symbol (due^, 
e.g., stands for the VEV of the color singlet combination of the components with one index having the 
quantum numbers of d, another of u and another of e. 
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(a) (b) 

Figure 2: Tree level diagram giving rise to the operator 0\ of Eq. fl3,9| ). All the indices should be considered 
as upper ones, except the ones for gauge bosons Q where upper and lower indices have been shown explicitly. 




Figure 3: Proton decay mediated by Higgs bosons, giving rise to the operator in Eq. ( 3.17| ). The notation 
about indices has been explained in Fig. One can similarly contemplate a diagram where the VEV of 
the component & ddUe appears instead of & due . 
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Figure 4: SU(15) symmetry breaking where baryon number is broken by the VEV of an antisymmetric 
rank-4 multiplet. The notation for VEVs has been explained in Fig. |]. 
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Figure 5: Tree level diagram giving rise to the operator O of Eq. (3.31). The notation about indices has 
been explained in Fig. |2[ One can similarly contemplate a diagram where the VEV of the component Q ddUe 
appears instead of <& due . 
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Figure 6: SU(15) symmetry breaking where baryon number is broken by the VEV of an antisymmetric 
rank-5 multiplet. The notation for VEVs has been explained in Fig. [I]. The two different VEVs of the 
multiplet H have been described in the text. 
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Figure 7: Higgs boson mediated diagram giving rise to the operator of Eq. (3.40). 
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Figure 8: Gauge boson mediated diagram giving rise to the operator of Eq. (|3,53| 
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Figure 9: Higgs boson mediated diagram giving rise to the operator of Eq. ( 3.56 ) 
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Figure 10: SU(16) symmetry breaking where baryon number is broken by the VEV of an antisymmetric 
rank-4 multiplet. The notation for VEVs has been explained in Fig. |]. 
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Figure 12: SU(16) symmetry breaking where baryon number is broken by the VEV of an antisymmetric 
rank-3 multiplet. The notation for VEVs has been explained in Fig. |l]. 
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